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1. INTRODUCTION 


Volume integrals associated with the integration of inhomogeneous 
Helmholtz equation are of practical interest in determining physical 
quantities in acoustic, electromagnetic and elastic fields. The inhomo- 
geneous scaler Helmholtz equation takes the following form: 

2 2 

7 4> + a $ = -4 ttp (Ij 

2 

where p(r) is the source distribution or density function, V and o are 
the Laplacian and wavenumber, respectively. It is well known [1,2] that 
a particular solution to Eq. (1] is: 

= jjj r'] R"^ exp iaR dV (2) 

n 

R = |r - r'l 

in which [4TrR)~^ exp icjiR is the free space Green's function, and n is the 
region where the source is distributed. The source distribution function 
p(^ can in general be either expanded or approximated in a polynomial 
form and hence P(^'5 is normally written as 

P(r') = (x'D^ (y')'^ Cz’)'’ (3) 

where A, v, v are integers. The integration of the vector Helmholtz 
equation is analogous, [2] . 

The reduction of time harmonic fields of frequency u in the acoustic 
and electromagnetic fields to that of integrating the inhomogeneous 
Helmholtz equation over a given volume can be found in many standard texts 
[3,4]. A formulation that leads to the required form of volume integration. 
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Eqs. C2.3) such that the elastic fields can be determined has recently been 
given in [5-7] . Using the dynamic version of the Betti -Rayleigh reciprocal 
theorem, an integral representation of the displacement field in an 
elastic medium containing an inhomogeneity can be given in terms of the 
eigenstrains 3-nd eigenforce ir? as: 


= - fff ^jkrs 

J J J 


TTj(r)dV 


(4) 


where g^^ are the spatial part of the free space Green's tensor function. 
For a linear isotropic elastic medium. 



1 r g2 ^ . exp i B R 

4t;p V R 

o 


- [■ 


exp ia R 


exp i 6 R , . 1 

-e-g— J.jn, I 


(S) 


in which is mass density, a and B are wavenumbers for longitudinal and shear 
waves, respectively. Expanding the eigenstrains and eigenforces in a poly- 
nomial of position vector r yields: 


n. (r) = A. + A.^ x^ + A.^^ x^x^^ + . . . 

= ®ij " ®ijk \ ••• 


(6a) 

(6b) 


and substituting it and (5) in (4), the displacement field is found to be 


u (r) = f . (r) A. + f .,,(r) A., + ... 
mx,-" mj v/ J mjk~-' jk 


+ F . .(r) B. . + F . .,,(r) B 
mi j So-* 13 mijkw 


ijk 


(7) 
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where 


^mj * 'f’mj ~ ‘^’mj 

2 2 

4ttp It) f -i (r1 = -6 - ~ ’i’l, 

o mjk'- k mj k,mj k,m] 


(8a) 

(8b) 


'^^Pq“ '"’m ^ij ^ ® '"’raij * 


^^Po“ = - [X a 5.. + 2u 8 

\,mij ^ 2u 


6 . 
mj 


(8d) 


Here, X, y are Lame's constants, and 
iff R ^ exp(ic 


i|;(r) 


exp(iaR) dV , 


(9a) 


^k(^^) = 


,-l 


x^R exp(iaR) dV , 


(9b) 


*^kl. . .s 


x^x^ ^ exp(iaR) dV, 


(9c) 


)(r) 


R“^ exp (i SR) dV, 


(9d) 


.(r) = 


,-l 


x^R exp (i SR) dV , 


(9e) 


^kl . . , s 


x^xj^ . . .XsR"^ exp(iSR) dV , 


( 9 £) 


This paper presents results for the volume integrals over a region 
that is either an ellipsoid, a finite cylinder or a rectangular parallelpipe 
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with semi-axes ai, a -2 and a^* Fig. 1. The integrals in (9) are subsequently 
referred to as the 4>-integrals and they are obtained in series form by expanding 
R ^ exp (iaR) in appropriate Taylor series expansions for regions r > r' and 
r < r', and by using the multinomial theorem with also the assistance of the 
classical result of Dyson [8] in the case of an ellipsoid. Certain derivatives 
of the ^-integrals that are of interest are also presented. 

Professors C. P. Yang and C. Saltzer of the Department of Physics and 
Mathematics at The Ohio State University, respectively, participated In many 
helpful discussions. 

2. SERIES REPRESENTATION OF THE ^-INTEGRALS 
Let R ^ exp ioR be expanded in a Taylor Series expansion for f' as 


R-' exp I.R 4 ^ [«: “P i»p] • ‘I") 

for r > r* 

and in a Taylor Series expansion for r as 

R-^ exp iaR =J^ ^ x. ^ [(r*)'^ exp iar'j . (11) 

for r < r' 

in which the summation convention is observed and i = 1,2,3. 

Employing the multinomial theorem as suggested in Ref. [1], 
the 4>-integrals can be explicitly written as triple sums: 


<» n n-il 

‘»>Cr) = I I I 

^ n=0 il=0 k=0 


(-1)^ 

V. K! (n-4-k)! 


.3 

. fc. k., n-Jl-k 
3x 5y 8z 




■ HI (»')' p(x’,y',z') dV’ , 

n 


for r > r* 


( 12 ) 
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and 


" n n /• 1 -.n , , , , 

. , , V V V C-1) i k Tl-l-k 

w = z z z fn-iirir r ' ' • * / ^ 

~ n=0 1=0 k=0 




3x' 3y' 3z' 
for r < r' 


( 13 ) 


The Taylor series representations given in Eq. (10) and Eq. (11) 
converge for the region r > r' and r < r', respectively. The 
integral (r) in Eq. (12) is normally used to evaluate physical quan- 
tities measured at large distance from the region n. The apparent singu- 

-1 -2 

larities present in Eq. (13) appear as iln e, e , e .... where e is a 
small positive number. These singularities disappear, however, if e is 
taken to be the radius of a sphere centered around the origin. In evalu- 
ating (r) for an ellipsoid, care must be taken in determining the con- 
tribution to the integral from the lower limit e. A further note on this 
is given at the end of Section 3. 


3. INTEGRATION OVER AN ELLIPSOIDAL REGION 


The integrals in (12, 13) are of either one of the following forms; 


4° = /// (x')P (y')^ (z')^ dV , 

n 


(14) 


4^ = /// P(x',y',z') 1 


|Sin^| 


dV 


(15) 


= III P(x-,y’,z-) , \ 


(16) 


5 



These integrals can be further evaluated as follows: 

= /// (x»)^ (y']® (z')^ dV 

n 

f+1 g+1 h+1 

^1 ^2 *3 4v RCf/2)R(g/2)R(h/2D 

(2m+ij C'2m+l) R^i) (17) 

0 If any one of the superscript power f,g, or h is odd, 
where a^^, a^, a^ are the axes of the ellipsoid, and 
2m = f + g + h 



R(m) 


(2m)! 

m! 


(18) 


This result was first obtained by Moschovidis [9] . 
(b) n = 0, 


<!> rr /// p(x',y', 2 ') 

n 


sin or' 


dV 


“ _ 1 2m-l _ _ 

= /// P(x’.y.z') I(-l) dV< 


00 2m-l 

_ Y , m- 1 ° e 

~ (2m -1)! \.p 

where 

\,p = (x*^+y^+z'V‘^ dV 


a 


X+1 u+1 v+1 
1 ^2 ^3 


(X+lj+v+2ni+l) 


4v 

(X+u+v+2m-l) 


(19) 
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I 

in^ f ^2 ^ ^3 


= 0 . 


2m, 2m_ 2m_ 

(m-1)! ^1 ^2 ^3 (2^l^^^X)! 

m^!ni 2 !ni 3 ! C,2m^+\)/2l ( 2 ni 2 +u)/ 2 I 

***1**"2***'3 ~ ’ 

if >.,ii, or V odd. 


(2ni2+w) ! (2nij+v) ! [2 Ctii.-1)+X+ii+v]/2! 
(2m^+v)72l [2(m-l)+X+y+v] I 

if X,y,v all even, (19a) 

(19b) 


in which the multinomial formula 

Cr')2"' = (x-2 ,y,2 ,,.2^r 


I „ m.r cy )"’"2 ( z .) 2">3 


"'l'"'2’”'3 3 


( 20 ) 


is used. In (20), the sums are taken over all non-negative integers m^. 


and mj for which mj^+m 2 +m 2 = m. 


(c) n = 0 , 4> 


^ = /// P(x',y',z') 

n 


cos or' 


dV 


(21) 


- I («•)' ty)" u')'’ ■ dv 


„ n (2m)! JJ 

m=0 ' 52 


Using the multinomial formula and letting the integral in (21) be denoted by 


^m,p " (r')^"*~^ dV 


V "■! Iff (x')2”'l^^Cy.)2"*2-Pf,,52m3+v 

^ m, !m_!m,! JT ■ dV» t J 

12 3 A 


m 


The volume integral in Eq. (22) may be viewed as the potential of variable 
densities observed at the origin, r=0. Applying the results on volume 
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integration over an ellipsoid given by Dyson [8], Eq. (22) can be 
written as 


m,p 


* 




(m)! 

™r ’"2- "s' 


2mj+X 2m2+v 2m3+v 



where 2p = (X+u+v) 



dx' 


2S.n 


dy'^^2 dz'^^'S 


( 24 ) 
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in which the sums are taken over non-negative values of for which 

P 

Zj^ + Z2 ■*' Z^ = Z. Using the definition of 6 , (24), and noting that 


, /a x'\^"’l'"^/a y> 

F= (x.)^"’l"\y)2"'2^^z.)2’"3"'’ \ /^2y ^ 




a2+'/< 


v: 

a2+^ 


2m3+v 


it can be easily shown that 


a"*"? F 


/v' V V’\ 

12 *2 ’ ~2 I 

\aj+4' a2+'l< 


(m^.X/2)!(m^^^p/lj !(m3W^ TT C2n>i-X) ! C2m2^p) ! (2m3.v) ! • 


1 I j \n.2W2/ j \">3W2 


aj+<|) 


2 , 

. aj+'l' , 




(25) 


Finally, for n = 0 


r 2m 

4>^ - y (-1)"' , C 

Jo ">>P 


(26) 


where 


TTEt! (2mj^+X) ! (2m^+y) ! ( 2 m 3 +v) ! 

“m,p ~ ^ „ m !m !m ! (m +X/2) ! (m_+u/2) ! (m_+v/2) ! 

m^,m2,m3, 12 3 1 2 3 
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2nn+X+l 2m2+p+l 2m3+v+l 

^2 ^3 


22m+2p (n,+p+i)i 


r dll/ 

J ^^ 2 ^^ jn . 2 + u /2 


( 27 ) 


(d) n ji 0, 


^ = ///p(x',y',2') 


sin ar' 


^ ,k„ ,n-Jl-k 

3x' 3y' 3z' 


dV 


00 

= I 

111=1 




2m-l 

(X 


(-1^ ('2111-1)! 



(28) 


where 


s" n = /// (x’)^(y’)''(z-)^ . 


m,p 


,n-£-k 
3x' 3y' 3y' 


. (r')2"'-2 dV' 




(m-D! (2mj)! (21112)! (2ni3) 1 
rm (n-!l-k) ! 


• III (x<)’^-^2mi-t ^y,^u+2m2-k ,^v+2iH3-n+i+k (29) 

n 

in Wiich the multinomial formula is used and m^, m2, m^ are summed over all 
integers greater than and equal to unity and mj^+m2+m2 are summed over all 

integers greater than and equal to unity and m^+m2+m2 = (m-1) . The integral 
in (29) can be obtained by using the formula given in (17), and is easily 
shown to be 
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47T 




-1)! (2mpi (2m2)! (2m^V. 

U k! Cn-A-k) ! 


u+2m2-k+l \>+2m,-n+l+k+l 

. ^1 ~ ^2 ^3 ^ 

(2p+2m+l-n) C2p+2m-n-l) 


R(5s(X+2mj-A+l))R(Jj(p+2m2-k-l))R(J2Cv+2mj-n+Jl+k+l)) 


R(Js(p+m-l-n/2)) 


if CX-il), (p-k), Cv-n+il+k) all even 


= 0 if (X-S,), (p-k) or Cv-n+Jl+k) is odd 


m,p 


I 

ra2«"'2’”*3 


(u 


X+2im -A+1 


where 


2p = X+p+v 


(e) n ^ 0, 


C rrr r , , I^ 3 cos or' ,,,, 

$ = ;/; pCx’.y'.z-) 1 ^ - n-a-k r’ -~ 

0, 3x* 3y’ 3z* 


2m 


I <-»"(liiTr 


m=0 


m,p 


where 




in.P .k^ ,n-Jl-k 

o 3x' 9y' 3z' 


Cr')^®’^ dV 


(30. a) 
(30. b) 


(31) 

(32) 
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When n 0, it is not as easy to find a compact form for these integrals. 

For the determination of the elastodynamic fields of an ellipsoidal inhomogeneity 

as formulated in [6,7] it is sufficient to determine $^(r) for a finite number 

of n's in determining the ^ , ®ijk’ ^jk’ example, 

♦ 

if it is necessary to determine the eigenstrains up to a second order 

(3 

distribution, it is then sufficient to find $ for 1 £ n £ 6. 

The integral $^(r) in (13) can be replaced for n = k by 


<)> (r) 

< .-V 


2 ! 


X X 

p q 


ao 


> 

m=0 


(- 1 )"' 


2m \ 
g p(k) 

(2m)! m,p 


(33) 


where 


m,p 


r f r 1 ^ 

= III (x') (y) (2') 

n 


8x' 3x’...8x' 
p q u 


dV 


(34) 


The substitutions of the derivatives of (r')^™~^ in Eq. (34), 
lead to integrals that can be easily evaluated by using 

Eqs. (23,24,25) , for the cases m > 1, k = 1, 3 and m > 2, k = 4 , etc. 

Special attention must be given to the cases m = 0 , k = 2,3, and m = 0,1 
k = 4. 

Using the notations given in Ref. [10] and noting that 

dV = dx' = dr' dS = dr'*r'^ dm , (35) 


we obtain 


r' W 



(36) 


where 



(37) 




g = l?/a? , 1. = x!/r’ 

and f = 0, e = 1 due to the fact that here we consider the source point 
being situated at the origin, i.e. _r = 0. Volume integrals associated with 
Eq. (32), 1 < n < 4, can be written into surface integrals by using Eq. (35). 

and finally reduced to simple integrals through the work of Routh [11]. e*g- 

///p*(r')“^ dV = f/p • (r')~^ • dr' • r’^ • du 

n 

~ III P ’ (r')~^ dr' dw 

If P = 1. [11, p. 901], 

/// (r')‘^ dV = 


/// Cr')'^ dV = 

n 


The surface integral of the type // ^3 reduced to simple 

integrals as well by using the work of Routh [7] in the same manner as 
listed in Ref. [6] and therefore will not be repeated here. The constants 


//{Hn r*(j..) - S,n e) dw 
S ^ 


- //{y iln g + in e } du 
I 


-l//{[r'Ci.)]‘^ - £-2} du, 

2 X 

Z 


(3 8) 


-1// g dm + 
^ Z 


2ir 


+ A- 


T a.a. "2 

X X 


(39) 
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_2 

The constants and are equal to 4Tr()l^a - and +(2ir/3)CE ), 

respectively for a sphere of radius a, where e is a small positive number. 

-1 -2 

The coefficient of these types of terms, in e, e , e , in the ^-integral 

I 

can be shown to be identically zero in a straight forward manner if n is a 
sphere. When n is an ellipsoid, the lower limit of integration should be taken 
from the surface of a small sphere with radius e, (38,39). The contribution 
to the 4-integral from the lower limit can therefore be identified as zero. 
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Fig. 1 An ellipsoidal region 
of integration. 
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